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We investigate the dynamics of colloids at a fluid interface driven by attractive capillary in- 
teractions. At submillimeter length scales, the capillary attraction is formally analogous to two- 
dimensional gravity. In particular it is a non-integrable interaction and it can be actually relevant 
for collective phenomena in spite of its weakness at the level of the pair potential. We introduce 
a mean-field model for the dynamical evolution of the particle number density at the interface. 
For generic values of the physical parameters the homogeneous distribution is found to be unstable 
against large-scale clustering driven by the capillary attraction. We also show that for the instability 
to be observable, the appropriate values for the relevant parameters (colloid radius, surface charge, 
external electric field, etc.) are experimentally well accessible. Our analysis contributes to current 
studies of the structure and dynamics of systems governed by long-ranged interactions and points 
towards their experimental realizations via colloidal suspensions. 
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I. INTRODUCTION 

In recent years the issue of structure formation by colloids at fluid interfaces has been the subject of intense 
experimental and theoretical research. If the colloidal particles are only partially wetted by the coexisting fluids, 
for not too high concentrations a two-dimensional (2D) colloid layer forms at the fluid-fluid interface because the 
detachment energy of a particle from such an interface is much larger than the thermal energy [l|, 0] ■ It turns out 
that these systems provide an excellent test-bed for fundamental issues (such as, e.g., 2D phase transitions 0]) as well 
as interesting perspectives for a variety of applications (such as, e.g., micropatterning and novel materials emerging 
from particle self-assembly Q). 

The colloidal particles forming a 2D layer at fluid interfaces interact with each other in diverse ways. For the 
present purpose, these interactions can be classified into three groups (see, e.g., Ref. Q)- (i) There is the ubiquitous 
van der Waals force, complemented by double-layer electrostatic interactions. These arc described by the Derjaguin- 
Landau-Verwey-Overbcek (DLVO) model (see, e.g., Ref. 0) and are relevant only if the particles have a chance to 
come sufficiently close, leading to coagulation, (ii) One can find also a repulsive interaction of longer range. The 
presence of unscreened charges on the particle surface exposed to a nonpolar fluid phase generates an unscreened 
dipole-dipole repulsion 041] ■ If both fluids are dielectric, the same kind of repulsion can be created by polarizing the 
particles with an external electric field [loj . If the particles are paramagnetic, a repulsion arises between magnetic 
moments induced by an external magnetic field Q. If one of the fluids is a nematic phase, director deformations 
induced by the anchoring boundary conditions at the particle surface can also lead to an effective repulsion [ll|, 113 ■ 
This second group of interactions is usually promoted in experiments with the aim to stabilize the colloids against 
coagulation, in which case one can effectively neglect the DLVO-type force, as we shall do in the following, (iii) 
Finally, upon integrating out the interfacial degrees of freedom one obtains an effective, capillary interaction due 
to the colloid-induced deformation of the fluid interface (see, e.g., Refs. [3, There are several ways how the 

particles can deform the interface. A capillary attraction arises if a vertical force is exerted on the particles, e.g., 
due to buoyancy p"5l - [l7j , an external electric field [l3| , or the effect of a substrate if the lower fluid is a film of finite 
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thickness (see, e.g., Ref. [iJl)- If the particles are non-spherical or the wetting properties of their surfaces are not 
homogeneous, the fluid interface is deformed anisotro pically and the ensuing capillary force is attractive or repulsive 
depending on the relative orientation of the particles |18l - [27| . 

The interest in capillary forces on the micrometer scale intensified recently when it was proposed [28| as an ex- 
planation for the puzzling attraction observed [l9l . [29l - [36| between micrometer-sized particles which were supposed 
to exhibit only the kind of repulsion discussed above as point (ii). However, various theoretical studies |37H4l| have 
shown that under the prevailing experimental conditions the capillary attraction between two spherical micrometer- 
sized particles is too weak to be able to explain the observed attraction. The reported effect might even turn out to 
be just an artifact of the experimental sample preparation [i^. Nevertheless, the capillary forces between charged 
particles still remains a topic of current research interest (see, e.g., Ref. [4^ for a recent study considering deviations 
from sphericity). 

These theoretical studies concern only the interaction between two isolated particles. The capillary interaction in 
the submillimeter range is formally analogous to unscreened 2D electrostatics (or 2D Newtonian gravity) [4l| and 
thus non-integrable in the sense of equilibrium statistical mechanics: the energy of configurations governed by this 
interaction is hyperextensive and the capillary attraction could be the relevant driving force for collective, genuine 
many-body phenomena in spite of its relative weakness at the two-body level. One of the goals of the present study is 
the investigation of this possibility for the kind of presently accessible experimental setups. Pergamenshchik [i^l has 
recently pointed out the possible relevance of this many-body effect as an explanation of the stability of the clusters 
observed in colloids at fluid interfaces. We shall show, however, that this claim is actually unfounded. 

The second goal of this work is to introduce a theoretical framework which addresses the dynamical aspects of 
collective evolution under capillary attraction. The corresponding studies published in the literature so far report 
on the motion of cither a single particle exposed to an externally created interfacial deformation or of few particles, 
usually two, following their own capillary attraction. Usually the equation of motion is solved in the overdampcd 
approximation in order to relate the velocity with the capillary forces, thus providing a means of interpretin g th e 
experimentally recorded particle position as a function of time (see, e.g., Ref. [45 1 and, more recently, Refs. [2^. I43ll46j ) . 
In Ref. (47j the full combined problem of particle motion and hydrodynamic flows in the fluid phases is addressed 
numerically with due account of the interfacial deformation. This approach is applied to clusters of two to four 
millimeter-sized particles. On that scale, the capillary attraction is screened and effectively very short ranged because 
it decays exponentially beyond the capillary length, which is typically of the order of millimeter. 

Beyond this interest in the detailed motion of individual particles it seems that little attention has been paid to the 
overall evolution of a colloid monolayer driven by its own, self-consistently determined capillary force field. In this 
respect we are only aware of Ref. (48| , where the clustering is studied experimentally as well as by means of molecular- 
dynamics simulations and interpreted within the framework of a certain effective kinetics of aggregation. In an attempt 
to be as realistic as possible, this latter theoretical approach takes into account many effects simultaneously (capillary 
forces, DLVO interactions, solvation forces, fluid streaming by temperature inhomogeneities) . As a consequence in this 
analysis the specific signature of the capillary interaction on the dynamics is masked. Here we study the evolution of 
the coarse-grained particle-density field, with the diffusion at the interface being driven by the interparticle capillary 
attraction at scales below the capillary length. This is addressed within the mean-field approximation as being 
appropriate for a non-integrable interaction at that range of scales. The problem has a formal resemblance with the 
evolution of a self-gravitating 2D fluid, which allows us to predict a phenomenology akin to the process of cosmological 
structure formation. 

The present study can also be viewed as a contribution to current investigations of the structure and dynamics 
of systems governed by long-ranged interactions, which can exhibit rather peculiar properties (iol. [50|. Our analysis 
points towards explicit experimental realizations of such systems via colloidal suspensions. 

In Scc.|TT]the mean-field model is introduced. After some qualitative considerations illustrating the overall picture, 
the model is formulated in terms of a set of coupled equations for the particle number density at the interface and 
the mean-field interfacial deformation. In Sec. IIIII the time evolution predicted by this model is analyzed within 
two different approximations which facilitate an analytical solution of the problem. We first investigate the linear 
stability of a homogeneous particle distribution and find a clustering instability, which is the analogue of the so-called 
Jeans' instability studied in the astrophysical literature. Then we compute the non-linear evolution of a radially 
symmetric perturbation in the limit of strong capillary attraction (the so-called cold-collapse approximation following 
the cosmological terminology). In Sec. IIVI we analyze available experimental setups and conclude that there is an 
accessible range of parameters for which one can expect that the predicted phenomenology is observable. In Sec.lVlwe 
discuss pertinent experimental observations reported in the literature in the light of our results as well as the precise 
relationship between our study and the work by Pergamenshchik [4^ . Sec. IVII summarizes our conclusions. 
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II. MEAN-FIELD APPROACH 



A. Qualitative considerations 




FIG. 1: (color online) Small deformation (exaggerated in this schematic drawing) of a fluid interface (solid line). M(r) is the 
vertical deformation field (i.e., in z direction) with respect to the reference, flat configuration (dashed line) as function of the 
lateral coordinate r = (x,y). n(r) is the vertical force per unit area acting on the interface in the reference configuration. 
These two fields are related by Eq. ([1]). 



We recall briefly the electrostatic analogy of small interfacial deformations, as worked out in Ref. [13, SH (see Fig.[T|). 
Let zi(r) denote the small vertical deformation of an otherwise planar interface at the lateral position r = {x,y), and 
n(r) the vertical force per mrit area exerted by external agents, e.g., a pressure imbalance across the interface or forces 
directly exerted on trapped particles (due to, e.g., buoyancy or optical tweezers). These two quantities are related by 
a Debye-Hiickel-type equation, 

T U 1 

V^u^-^--U, (1) 

where 7 is the surface tension and A is the capillary length. This equation describes local mechanical equilibrium: at 
each point of the interface, the force by external agents (H) is counterbalanced by the surface tension of the curved 
interface (V^u), augmented by the force due to the weight of the fluid displaced relative to the flat conflguration 
(it/A^). The lateral force on a piece S of the interface exerted at its border by the rest of the interface takes the 
simple form 

Flat - f d^rIl{r)Vu{r). (2) 
Js 

Equations ([1]) and ([2]) lead to the analogy that the deformation u plays the role of a 2D electrostatic potential, H is 
the charge density, and the capillary length^ A is the screening or Debye length, with the peculiarity that the forces 
are reversed and charges of equal (different) sign attract (repel) each other. Or, in other terms, the analogy holds 
with a "screened" gravitational interaction involving positive and "negative" masses. 



The capillary length is given by A = \/^/{g\pi — P2 1) in terms of the acceleration g of gravity and the mass densities p\^2 of the coexisting 
fluid phases. 
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FIG. 2: (color online) Side view of a fluid interface containing many particles, each subjected to a vertical external force /. 
Within the simplest approximation the particles are modelled as pointlike capillary monopoles of strength /. 



In this respect, a particle trapped at the interface is characterized by a set of (capillary) multipolar charges^ and the 
lateral force is expressed in terms of the coupling of these multipoles with the deformation field u{r). Thus, similarly, 
the capillary interaction between a collection of particles (sec Fig. [2]) can be written as the sum of the interaction 
between pairs of multipoles localized at the position of the particles. In particular, the (isotropic) capillary monopole 
is simply equal to the net vertical force / exerted on the particle by external agents, i.e., aside from the force provided 
by the interface [4l|, [s^] . The effective interaction between two monopoles of equal strength / separated by a lateral 
distance d < A is described by the attractive potential [l^-IT7| 



V{d) « --^ 



27r7 a 



(d < A). 



(3) 



(Actually, the exact potential is given by V{d) = —pKo{d/X)/{2Tr'-f) in terms of a modified Bessel function; accord- 
ingly the length scale in the logarithm in Eq. ([3]) is w 1.12A instead of A, but we neglect this small difference for the 
qualitative estimates to follow. At separations d ^ A, the interaction is screened and the potential crosses over to an 
exponential decay, which, for reasons of simplicity, we set to zero in the following reasoning.) This a 2D Coulombic 
interaction, which is known to be non-integrable and to drive an instability, e.g., in the astrophysically relevant case 
of a self-gravitating gas. 

At this point we consider it to be useful to introduce a qualitative discussion which illuminates the physical origin 
of this instability. The more formally interested reader can skip this part without loss in favor of the mathematical 
derivation presented in the following sections. Consider a collection of iV ^ 1 identical monopoles at positions 
(i ~ 1, . . . ,N) distributed homogeneously over a region of linear extension L (and thus of average number density 
g ~ N / L"^). Due to the long-ranged nature of the interaction, the capillary energy per particle Ccap of the configuration 
can be estimated at the mean-field level as follows: 



—yv{\i 

2N 4^ ^' 
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(4) 



The key point is that for L < A this energy is not an intensive quantity but scales instead like N = gL^ . There is also 
a contribution eghort to the energy per particle, due to thermal motion and short-ranged, predominately repulsive 
forces, which is TV-independent. Therefore, if the system is large enough (but still L < A) it can happen that the 
absolute value of the energy due to capillary effects Ccap (< 0) dominates over eshort (> 0). The corresponding critical 
system size can be characterized by Jeans' length Lj (this terminology is borrowed from the astrophysical literature), 
which is determined roughly by the condition Cshort — le^ 



^cap I 



1 / 87eshort 



/V g 

upon neglecting the logarithmic correction in Eq. (|4]). (A precise definition will be given in Sec lIIIi see, c.f., Eq. ([19 



(5) 



^ If R denotes the characteristic size of the particle, the multipolar expansion corresponds to, strictly speaking, an intermediate asymptotics 
at distances R <^ r ^ X; for r ^ X, the expansion can be re-summed (see, e.g., Ref. ISlll ). The corrections to the multipolar expansion 
are, however, suppressed by powers of the small ratio R/X. Since R ~ 1-10 /im and A ~ 1mm in the problems we are interested in, we 
can safely neglect this effect. 
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The capillary energy (Eq. (|4])) can be expressed in terms of this length as follows: 

L<X, 
\<L. 



X 

Ccap — ~^shoit I "7 I ^ 











1 - 


h21n^' 






A 



(6) 



One can distinguish two distinct extremal cases (see Fig. 

• A < Lj: One can easily see that in this case |ecap| < Cghort, independently of the lateral extension L of the 
homogeneous particle distribution and the effect of the capillary attraction is just a small perturbation. This 
case corresponds to the physical situation in which the capillary force is screened (i.e., negligible at distances 
beyond the capillary length A) preempting that its cumulative effect becomes comparable with the net effect of 
the short-range forces associated with eghort- 

• Lj < X: Here, |ecap| > Cghort whenever Lj < L, implying that the homogeneous distribution is unstable because 
the attractive capillary force cannot be counterbalanced by the nonzero compressibility provided by Cshort- The 
system collapses into an inhomogeneous, clustered state with a new typical size L* and a new average density 
g* sustained by its own capillary attraction energy e*ap, and for which the compressibility provided by 6*)^^^.^ 
can balance the capillary attraction so that |e*j^p| ~ e*j^Qj,^ (and the new, effective Jeans' length Lj is given by 
Eq. (O). If Eq. ^ is evaluated for the quantities carrying an asterisk (characterizing the clustered state), this 
condition has two possible solutions: 

(i) the cluster size is of the order of its effective Jeans' length, and the latter in turn is smaller than the 
capillary length, i.e., L* ^ Lj < A, 

(ii) the effective Jeans' length is of the order of the capillary length, and the latter in turn is smaller than the 
cluster size, i.e., L*j ^ X < L* . 

(In Appendix [C] we present an alternative qualitative derivation of condition (i) based on a force-balance 
argument rather than this energy consideration.) 

These qualitative considerations are formalized in the following sections using a mean-field model for the particle 
dynamics under the action of capillary attraction. 





FIG. 3: (color online) Qualitative stability diagram of the homogeneous distribution in terms of the capillary length A and 
Jeans' length L,j as a function of the system size L for the two cases X < Lj (a) and Lj < A (b). 



B. Mean— field capillary force 



We consider a system of particles trapped at a fluid interface considered to be flat for the time being. In the absence 
of long-ranged interactions, there are equilibrium states corresponding to a 2D homogeneous fluid phase with areal 
particle number density characterized by an equation of state p{g), where the "pressure" p is the lateral force per 
unit length exerted by the particles on the walls providing the lateral confinement of the monolayer. (In what follows, 
we omit the possible dependence on the temperature T which is irrelevant in the present context because we assume 
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isothermal conditions maintained by the fluids on both sides of the interface.) One can write (setting the Boltzmann 
constant equal to unity) 

p{q) = qT +Pc^{q), (7) 

where the first term is the entropic (ideal gas) contribution and the second term (excess pressure) is determined 
by the short-ranged interactions between the particles. There are always the hard-sphere contribution and the van 
der Waals attraction which eventually leads to colloid coagulation. However, in the experiments aimed at observing 
capillary effects these contributions are irrelevant because a longer-ranged repulsion (of electric, magnetic, or elastic 
origin) is implemented precisely to avoid coagulation [sj . Typically this repulsive potential has the simple form 

^rep(d) = = T (^0 , {n>2), (8) 

where the constant potential parameter A can be replaced by the Bjerrum length C, = (A/T)^/" of the potential (i.e., 
v,-c-p{d = 0= T ^ so that Q decreases with increasing T). The repulsion of electric or magnetic origin corresponds to 
n = 3 @, 3 : while the value of C, can be varied via the externally controllable parameters of the system (strength of 
the external magnetic or electric field, water salinity, temperature, etc.). Because of the simple scaling behavior of 
this potential, the phase diagrams of such kind of fluids are determined by the single dimensionless parameter C,"^ g 



(see, e.g., Ref. l53|): for sufficiently small C^g (high temperature, i.e., small C,, or low density) there is a fluid phase; 
if C^g is above a certain value, depending on n, the system freezes (for n = 3 this threshold value is C^p w 4.6 as 
obtained from experimental data [3|). Since the system is two-dimensional, the fluid-solid transition is expected to 
be of the Kostcrlitz-Thouless type. For n = 3 this is supported experimentally 0] and by recent simulations 54 1. 
according to which there is a narrow range of values of the parameter C,'^ g within which a hexatic phase is observed 
between the fluid and the solid phase. 

In order to address the effect of the long-ranged capillary force, one introduces the ensemble-averaged interfacial 
deformation U{y). From Eq. ([T|) one obtains 



V^C/-^ = -4, (9) 



in terms of the particle number density field g{v). Here, the ensemble-averaged vertical pressure H(r) has been 
replaced by the density of capillary monopoles, 

H(r) = H^o„(r) - fg{v), (10) 

where the capillary monopole / associated with a single particle is the net vertical force exerted on it by external 
agents. By analogy with 2D electrostatics and gravity, one can apply a mean-field approximation: the lateral capillary 
force experienced by a single particle located at r is written as +/VJ7(r), after replacing the ensemble-average of 
nVu in Eq. ^ by fgWU. Within this approximation correlations are neglected because U{r) is computed from the 
field g{r) via Eq. rather than from the density field conditional to the presence of a particle at r. Therefore U{r) 
is actually the coarse-grained correlate of the interfacial deformation u(r), neglecting small-scale spatial variations. 
The physical assumption underlying the mean-field approximation is that the dynamics of a particle is predominantly 
determined by the simultaneous interaction with many other particles. This is usually expressed in terms of the 
constraint that the parameter gX^ number of neighbors which de facto exert a force on a particle) must be large. 
In the experiments of interest here this is always the case, because A ~ 1 mm but the mean interparticle separation 
(~ g^^^^) lies in the micrometer range. (Another implicit assumption, peculiar to the capillary problem, is that the 
capillary monopole / of a particle is independent of the presence of other particles, i.e., of the particle density g. 
This is a good approximation if the vertical force / is predominantly due to gravity or due to an external electric or 
magnetic field, see, e.g., Ref. Q and references therein.) 

In summary, the equilibrium state of the system is described macroscopically by the force balance equation 

- Vp(^?) + /^.vc/ = 0, (ii) 

which, together with Eq. ([7]) and Eq. determines the equilibrium density profile g{r). (Note that Eq. ([TT]) also 
follows from, c.f., Eq. ([H)) . or more generally Eq. (jASp . for dg/dt = 0.) Except for the finite value of the parameter A 
in Eq. this problem is formally analogous to the Vlasov-Poisson model (corresponding to A = oo) for determining 
equilibrium configurations of a fluid under its own self-gravity [ssf (see Ref. (56j for a comprehensive analysis of the 
equilibrium solution of the Vlasov-Poisson model for arbitrary spatial dimensions with Pexiff) = 0, i.e., for an ideal 
gas, see Eq. ([7])). In recent years this old problem has received renewed attention and generalizations of it have been 
studied throughly (see Ref. [s^l and references therein for a brief summary). 
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C. Diffusive dynamics 

If the force balance equation pT|) is violated, the particle density p(r, t) will evolve in time according to the law of 
mass conservation, expressed by the continuity equation 

|=-V.(.v). (12) 

The flow velocity field v(r, t) is determined by the law of motion of the particles. Each particle is dragged at the fluid 
interface by a mean force — (Vp)/p + /VJ7. We assume that the characteristic time scale of macroscopic evolution is 
long enough so that the motion of the particles at the interface occurs within the overdampcd regime. This allows 
one to neglect particle inertia and a flow velocity is induced given by 

v = r(^-^ + /vc/^, (13) 

where F is a mobility coefficient of the particles at the interface. In addition there are hydrodynamic interactions 
between the particles due to the fluid flow induced by the particle motion which can affect the dynamical evolution 
[ssf . The effect of this interaction could be incorporated through a density dependence of T (see, e.g., Ref. [s^ for 
sedimenting hard spheres in bulk fluids) or, in a more explicit manner, by additional terms in the diffusion equation 
p4p below [60j . For our purposes, however, we neglect this effect and consider F as a phenomenological input 
parameter which is taken to be spatially constant for reasons of simplicity"^. Inserting this flow field into Eq. (I12p one 
obtains 

^=rW-[Wp{g)-fgWU]. (14) 

On the other hand, we assume that the evolution of the areal number density profile g{r, t) occurs on a time scale 
sufficiently large so that deviations from local equilibrium and the presence of capillary waves can be neglected. 
Therefore, the equilibrium relationships in Eqs. ([7]) and ^ hold and, when combined with Eq. (jl4p . a closed equation 
is obtained to determine the shape and the evolution of the density distribution g{v, t) , which thus follows a diffusive 
dynamics driven by "self-gravity" . As discussed in Appendix |^ the problem can be cast in terms of a functional 
formulation. 



III. CLUSTERING INSTABILITY 



A. Linear stability of a homogeneous state 



The occurrence of a clustering instability can be inferred from Eq. (|14p . For a macroscopically extended interface 
a homogeneous particle distribution of density gu is a solution of Eqs. ([9|) and dTTI) with Uh = f^^Qh/l ( = const.). 
We consider now a perturbed configuration g(v,t) = gn + 5g[r,t), U{r,t) ~ Uh + SU{r,t) and linearize Eqs. © and 
(fH)) in terms of the small perturbations Sg and SU: 



djSg) 
dt 



FV 



QhK-h 



-Wdg^fghS76U 



(15) 



6U 



-Sg, 



(16) 



where the isothermal compressibility Kh is given by 



dp 



l^h ■= ( g-g^j {Q = Qh) 



(17) 



^ This approximation holds in the dilute limit. Two particles of size R moving in a three-dimensional bulk fluid at a distance d ^ R 
acquire a relative velocity correction of the order (R/d) |61[ . On the other hand for a solution of sedimenting hard spheres of radius -R 
and bulk number density gi,, the mobility F is corrected by a factor fa 1 — (3R)^ Qi, to lowest order in gj, ISSH . These results can serve 
as a first estimate of the effect of the hydrodynamic interactions. However, one should keep in mind that the computation of these 
interactions in the presence of a deformable interface is actually an open problem, which lies beyond the scope of the present study. 
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We introduce the spatial Fourier decomposition of the perturbation, 



(18) 



and define a characteristic wavenumber K and a characteristic time T associated with the unperturbed homogeneous 
distribution as 



Polish 



T: 
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The diffusion equation reduces to 



dSg 

dt 



r(fc) 



with a typical time of evolution T(fc) given by 



T(fc) 



QhUh 7 1 + (Afc)^ 



1 //c 



t/r(fc) 



{k/KY + (AA')- 



(19) 



(20) 



(21) 



In the astrophysical literature (in which A = oo for gravity, see, e.g., Refs. |55l.l62| and references therein), K is known 
as Jeans' wavenumber (and Lj :— 1/K is the associated Jeans' length, see Eq. ([5])). The value of this parameter 
is determined by the properties of the unperturbed homogeneous state. Two qualitatively different cases can be 
distinguished (see Fig. S]): 

• \K < 1, so that T{k) < for all values of k. In this case perturbations of all wavelengths decay exponentially 
as function of time; therefore the homogeneous state is stable. Physically, this describes the situation in which 
the number of particles inside the circle of interaction of radius A is too small and thus the capillary attraction 
is too weak to lead to a collapse of the colloidal fluid against a finite compressibility. 



AA' > 1, so that T(fc) > for wavcnumbers below a critical one. 



Kjl- 



[XRy 



(22) 



determined by the condition T(kc) ~ 0. Perturbations satisfying this condition arc linearly unstable, which 
describes the onset of a clustering instability. In the limit of no screening of the capillary attraction (i.e., for 
A — >■ oo), one recovers the scenario of Jeans' instability: any homogeneous state is unstable against perturbations 
with a wavenumber smaller than Jeans' wavenumber. Figure [5] depicts kc as function of A: for all practical 
purposes one can take kc ~ K unless the parameter AA' is close to one. As can be inferred from Fig. 21 T is the 
characteristic time of the fastest growing mode if AAT is not too close to one. 

The value of Jeans' length Lj ~ 1/K associated with every homogeneous configuration determines its stability 
against clustering by capillary attraction. If L denotes the linear extension of the system, the results are summarized 
by Fig. [3] obtained previously in Subsec. Ill Al based on qualitative arguments. The dependence of stability on the 
equation of state of each particular system enters only through the definition of Jeans' length (Eq. P^ ). A clustered 
phase is only possible if Jeans' length is small enough, AAT > 1, and the system size is sufficiently large, i.e., 27r < LK, 
because a system of linear extension L cannot support perturbations with wavenumbers smaller than 27r/L. (Actually, 
if the system has a finite extension, the theoretical analysis above has to be complemented with appropriate outer 
boundary conditions compatible with a homogeneous particle distribution. We do not expect this to change the 
conclusions qualitatively, but maybe the precise value of AAT setting the boundary of the stable region in Fig. [3l) 

Finally, we remark that this Jeans-like instability has also been analyzed in two recent studies. In Ref. 16311 a model 
of bacterial chemotaxis is considered, to which our equations reduce formally in a certain limit. In Ref. |64| the same 
mathematical model as ours is studied except that the system is confined to a disk-shaped region with Neumann 
boundary conditions. 
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FIG. 4: (color online) Inverse relaxation time l/rik) (Eq. as function of the wavenumber of the perturbation for different 

values of the parameter \K. The units T and K are given by Eq. (|19p . For \K < 1 one has T{k) < and thus linear 
stability against perturbations of all wavenumbers. For XK > 1 perturbations with long wavelengths are linearly unstable due 
to r(fc) > 0. 



Ir 




FIG. 5: (color online) Dependence of the critical wavenumber kc (Eq. ((22])) on the capillary length A. 



B. Cold collapse of a radially symmetric perturbation 

Going beyond this simple linear analysis is possible only by resorting to numerical computations. There is, however, 
a case which can be addressed analytically: the dynamical evolution of a radially symmetric perturbation of a 
homogeneous background in the limit of vanishing Jeans' length, meaning physically an arbitrarily compressible fluid 
(sec Eq. (jl9p ). In the cosmological literature this scenario is termed a "cold collapse" because in this context one 
considers ideal gases; for them an infinite compressibility amounts to a vanishing pressure p{g), which corresponds 
to the limit of zero temperature. The cold collapse is therefore the limiting case of a more general scenario involving 
both the capillary (gravitational) attraction and the pressure opposing compression. This approximation allows one 
to obtain an exact analytical solution of Eqs. ([9|) and (jl4[) in the presence of radial symmetry in the limit A — >■ oo. The 
computational details are presented in Appendix [C] Here we just summarize the results: The evolution of a localized 
radially symmetric perturbation of a homogeneous configuration is driven predominantly by the capillary attraction 
provided the system size is larger than Jeans' length; in this case collapse occurs only if the spatially averaged density 
of the perturbation is larger than the one of the homogeneous start configuration (i.e., if there is an overdcnsity) . 
The collapse proceeds until a cluster size of the order of the cluster Jeans' length (i.e., the one associated with the 
density of the cluster) is reached and the pressure is able to halt the collapse. This last stage of the evolution, which 
does involve the effect of a finite compressibility, is beyond the cold-collapse approximation by construction and we 
have not studied it yet. The total time of collapse is given by Eq. (|C13I) within the cold-collapse approximation and 
is roughly of the order of the characteristic time T associated with the homogeneous configuration (see Eq. ([T^). 

This exact solution of a simplified model can be used to gain some insight into a more reaUstic situation. Suppose 
that in a homogeneous configuration of density Jeans' length is small enough so that there can be found patches 
of linear extension L < X fulfilling the condition of instability, Lj < L (see Fig. [3]). In such a patch there are 
always thermally induced density fluctuations dg. Far from phase transitions these exhibit a Gaussian distribution 
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and are uncorrelated. If such a patch is large enough to contain many particles, the general theory of thermodynamic 
fluctuations states that the relative amplitude of these fluctuations scales like the inverse of the total number of 
particles in the patch, i.e., \5g\/ Qh — '^1 ^ 1- The cold-collapse model can be employed to estimate the time 

of collapse in that region with such an increased density (see Eq. (jC13p with go = Qh + l^oD- 

TcouW^riny^. (23) 

Since Tcon{L) increases with L, one would expect a bottom-up scenario of cluster formation in the language of 
cosmology, according to which spatially smaller perturbations collapse first, as opposed to a top-down scenario. 
However, the weak logarithmic dependence implies that the bottom-up clustering would be hardly observable; it is 
likely to observe the almost simultaneous collapse of fluctuations of all sizes into a single cluster of maximal size. 



IV. FEASIBILITY OF EXPERIMENTAL REALIZATIONS 



In the previous sections we have shown that the instability is characterized by two parameters: Jeans' length 
1/K and Jeans' time T (see Eq. ([T^ V In this section we compute these parameters for different setups which can 
be realized experimentally. The clustering instability will be easily observable if one can find a range of physical 
parameters for which the following constraints hold simultaneously: (i) the particle size R and Jeans' length 1/K 
should satisfy R,l/K < A (Fig. [3]), (ii) the mean interparticle separation q, which in this work will be measured in 
units of the particle size R, i.e., 

9:=7^, (24) 

should be small enough so that X/{qR) ^ 1 and the mean-field predictions apply, and (iii) Jeans' time should lie 
within a reasonable range which permits observations following up the collapse. At this point we emphasize that in 
the present context the goal is to observe a collective effect, i.e., a many-particle instability: Capillary attraction is 
routinely observed between particles visible for the naked eye (i? > 1 mm) . but in such a case the capillary attraction 
is actually a force of very limited range (^ A) and the corresponding phenomenology is completely different. Here, 
however, we focus our attention on micrometer-sized particles so that i? <C A. 

The values of Jeans' length and Jeans' time depend on the specific physical system via the strength / of the capillary 
monopole and the compressibility k determined by the equation of state p{g) (see Eq. P^ ). We shall consider in 
detail three different systems which are customarily employed in experiments: (i) Charged particles or particles with 
dissociable surface groups, for which / is due to their weight and p(g) is determined by the electrostatic interaction 
between the particles, (ii) Neutral particles at the interface between two dielectric fluids in the presence of an external 
electric field. In this case the capillary monopole is due to both the weight of the particle and the electric force exerted 
by the external field, while p{g) is again determined by the electrostatic interparticle force, (iii) Superparamagnetic 
particles in an external magnetic field. Here the capillary monopole is due to their weight and the 2D equation of state 
is determined by the magnetic interaction between the particles. In the following two subsections, we first compute 
Jeans' time, which depends only on the monopole /, and then Jeans' length, which in addition involves the equation 
of state. 

In order to be specific, in the following we consider typical values 7 = 0.07 N/m (surface tension of the air-water 
interface at room temperature) and l/F = GnrjeffR, where the effective viscosity rjes interpolates between the viscosities 
of the adjacent fluid phases (see, e.g.. Fig. 4 in Ref. [1^ and the measurements reported in Ref. [i^) and depends 
on the contact angle at the particle-interface contact line; we take 77eff = ?7water/2 = 0.5 x 10"'^ Pa x s (corresponding 
to a sphere half immersed in water at an air- water interface at room temperature). Because of the simple scaling of 
Jeans' length and time with 7 and F (see Eq. ([TO)) ) one can easily obtain the values of 1/K and T for other values of 
7 and F from the estimates we shall quote below. 



A. Jeans' time 



1. Monopole due to buoyancy 



We first compute Jeans' time T (see Eq. ([T^ V This time depends only on the strength / of the capillary monopole 
and is independent of the detailed form of the interparticle repulsion. Every particle at the interface carries a capillary 
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monopole / due to its weight (corrected for the buoyancy effect due to the fluids) . A spherical particle of radius R 
floating at a fluid interface experiences the vertical force 

47r _ 3 

/buoy = -^gPoff-R , (25) 

where g is the acceleration of gravity and poff is a (signed) effective mass density, which depends on the mass densities 
of the particle and of the fluids. For an estimate, we take as a typical value pcS ~ — 1-6 g x cm^^ (corresponding to 
the glass particles at the interface between air and corn oil used in the experiment described in Ref. [6^), so that 

Auo.«-i^f^)', (26) 
/im \ /im / 

with Troom = 300i^ as the room temperature and the minus sign indicating that the force points downwards. Fig- 
ure inja) shows Jeans' time as a function of the radius R for several values of the average interparticlc separation q. 
In this case, Eq. yields the simple scaling T cx q^R^^. 



2. Monopole due to external electric field 



Alternatively, a capillary monopole can be generated by a vertical external electric field E which polarizes the 
particles at the interface between two dielectric fluids [l3,[6^^. The vertical force on a spherical particle of radius R 
due to such a field E is 

/eicc = £oi?'£^'(e2-ei)0, (27) 

where ei and 62 are the dielectric constants of the upper and lower fluid, respectively (if the electric field points 
upwards) and the factor (j) depends on the dielectric constants and the height of the particle at the interface (see 
Fig. 7 in Ref. [l^, where the factor^ is called /„). With the values ei = 1, ea « 2.87, (j) « 0.27, and E 10^ V/m for 
the experiment described in Ref. [66[, one obtains 



/im \10^V/m/im 



and the total capillary monopole is the sum 

/ — /buoy 4~ /clcc ^ 



-^room / R 



Hm \/Ltm 



10^ 



16R 



lO^V/m/ /jm 



(29) 



Figure [6l^b) shows Jeans' time as a function of the radius R for several values of the electric field for a fixed particle 
density corresponding to a mean interparticle separation q = 10 (in units of R). In this case, the electric field 
can compensate the weight and a neutral buoyancy (i.e., / = 0) can be achieved at a specific value i?*//im w 
[E/ (1.26x 10^ V/m)]^ of the radius for a given electric field. Close to this value, Jeans' time diverges as T ~ (i?*— i?)~^. 
For radii much larger than this critical value, T approximates the buoyancy dominated regime discussed before, while 
in the limit R <^ R^, the capillary monopole is dominated by the electric force, i.e., / w /dec and Toe q^R~^E~^. 



3. Monopole due to external magnetic field 



In an experimental setup similar to the one just considered, superparamagnetic particles can be placed in an external 
magnetic field H perpendicular to the interface which induces a capillary monopole due to the magnetic vertical force 
which is described analogously to Eq. ([27| . However, the small values of the magnetic susceptibilities xi and X2 



The dipole induced in a particle creates an electric field E^ip{r) decaying far from a particle asymptotically ~ r~^. The action of 
this field on the interface induces an additional deformation not addressed in these studies. It can be computed from Eq. |[T} with the 
electric pressure 11 oc (i? + E^ip(r))^ , leading to an additional intcrfacial deformation u r~^. Here we neglect this contribution in the 
asymptotic comparison with the monopolar deformation, but acknowledge its potential relevance at high particle densities. 
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of the upper and lower fluid phase, respectively (typically |xi,2| ~ 10 ^) renders this force irrelevant under usual 
experimental conditions. Similar to Eq. (|27p one can obtain the estimate (up to a geometrical factor of order 1) 

„«^H»(^-^).i2:!Z=(^ A) ^ ,30, 

V1 + X2 1 + Xi/ Aim \ 10"^ A/m ^m / 

which, even for strong magnetic fields, is much smaller than the capillary monopolc due to weight given by Eq. (|26l) . 
Therefore, in the following we shall no longer consider the effect of a magnetically induced monopole^. This notwith- 
standing, in the following subsection we shall study the influence of an external magnetic field on the compressibility 
of a monolayer of superparamagnetic particles and thus on Jeans' length for gravity induced capillary monopoles. 




FIG. 6: (color online) Jeans' time for a homogeneous distribution of spherical particles as function of its radius R. The particle 
mobility at the interface is taken to be (see main text) F = {3ttR)~^ m x s/g and the interfacial surface tension is chosen as 
7 = 0.07 N/m. The 2D number density qh = l/(qR)^ is parametrized in terms of g, i.e., the mean interparticle separation 
divided by R, see Eq. (|24|l . In (a) the capillary monopole is given by Eq. (|26p . i.e., under the sole action of gravity, and 
corresponds to a vertical buoyancy force directed downwards. In (b) the capillary monopole is given by Eq. (|29fl . i.e., the 
particles are under the combined action of gravity and of an external electric field. The mean interparticle separation in units 
of R is fixed at a value g = 30 and the reduced electric field E is -5/(10® V/m). The divergence oi T at R ~ 15.7 /xm for 
E = 0.5 corresponds to a vanishing total capillary monopole. For E — 1 and E = 2 this divergence occurs at -R ~ 63 ^m and 
R ~ 252 /im, respectively. For this reason the curve for E = 1 exhibits a minimum on the scale of the figure. 



B. Jeans' length 

1, Compressibility for repulsive interactions 

The determination of Jeans' length requires the specification of the equation of state of the 2D fluid of colloidal 
particles (see Eqs. ([7]) and ([T9| ). The simplest case is that of an ideal gas, Pox(e) = 0; the corresponding Jeans' length 
is given by 



1 



^ ideal 




(31) 



Here the capillary monopole is only due to buoyancy so that / = /buoy (see Eq. ([26])) and K^^-^ cx qR^^. This length 
is plotted in, c.f.. Fig. [Slja) for reference. Similarly, one can consider the reference case of a 2D gas of hard disks of 
radius R (which coincides with the radius of the spherical colloid al p articles if they are half-immersed in one of the 
fluids). The equation of state is described well by the expression [63 

p^gT^^, (32) 



^ In the experiments described in Refs. [3I. ISSH . the particles are completely wetted by water and thus remain submerged, but very close 
to the interface. However, it is conceivable that superparamagnetic particles can be prepared which arc only partially wetted by the 
fluid phases and thus get trapped at the interface. 



13 



where gc := l/{2\/3R^) is the number density for close packing of disks in 2D. The corresponding Jeans' length can 
be written as 

' - (33) 



K Ki. 



ideal 



where 1/A'idoai is Jeans' length of an ideal gas at the same temperature and with the same number density (see 
Eq. (131])) and 




is a dimensionless function with £(0) = 1 which collects the deviations from the ideal gas behavior. After taking 
Eq. ([5^ into account this correction of the ideal gas behavior is a function C{R^g = q"^) solely of the dimensionless 
parameter q defined in Eq. (j24p . which must be larger than its value Qc — ^/{R\/q2) ~ 1.861 at close packing. One 
finds that this correction is significant actually only very close to q = qc] e.g., for q ~ 3 one has C ~ 2.07 so that in 
practical terms the curves in, c.f., Fig. [8l^a) are applicable also for hard disks. 

However, neither an ideal gas nor a gas of hard disks correspond to the generic experimental situation. Typically 
the particles are endowed with a soft interparticle repulsion in order to avoid coagulation brought about by attractive 
dispersion forces. This kind of repulsion is described by Eq. (|8]), where n = 3 corresponds to the experimentally 
relevant situations we want to address. As stated in the context of Eq. ([5]), for such interaction potentials the phase 
diagram depends on the parameter C^p only. We have used MC simulations (details can be found in Appendix |B]) to 
compute the e quat ion of state of this specific 2D fluid (see Fig. [7{a)), for which we are not aware of published data 
(compare Ref. [53] for other values of n in a 3D fluid). This equation of state is valid only for <^^q ^ 4.6; beyond 
this value the 2D liquid freezes into a solid phase Jeans' length can be expressed again as in Eq. but the 

function C defined by Eq. ([M)) is now a dimensionless function of the parameter ("^g (see Fig.[71Jb)). At low densities 
there is a relatively weak divergence, 1/K oc l/y^ (see Eq. (PTj) '). reflecting the increase of Jeans' length caused by 
a weakening of the overall capillary attraction upon dilution. At high densities close to the onset of freezing Jeans' 
length increases slowly: l/K - g^/^ because £ - (C^e)^/^ as provided by a numerical fit and due to l/i^idcai ^ 1/%/^ 
(see Eq. (ISg)- 




FIG. 7: (color online) 2D fluid of particles interacting according to the pair potential given by Eq. ((8]) with n — 3. (a) Pressure 
p as a function of (^'^g, so that p/{gT) — > 1 for p — s> 0. The dots are the results from a MC simulation, while the lines correspond 
to the approximation of the equation of state up to the second term in the virial expansion, p « gT[l + {■k/2)T{1/3)(^'^ g], valid 
for low densities, and to the numerical fit p ~ gT[3 + 6.6{^^ g)^^^], valid for large densities. Note that this latter fitting function 
corresponds to the equation of state of an harmonic solid of soft particles interacting with the pair potential [g^. We 
have found that in good approximation this equation of state also holds for the corresponding high-density fluid phase. For 
> 4.6 the fluid fr eezes. (b) The dimensionless function £ given by Eq. p4|) which characterizes Jeans' length (Eq. (|33[l V 



2. Repulsion between electrically charged particles 



A common way of implementing the interparticle repulsion is to use electrically charged particles or, more frequently, 
cover the particle surfaces with chemical groups which dissociate in water. (Even then, the repulsion is long-ranged 
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if the adjacent fluid phase is a dielectric 0.) The repulsion is described by Eq. ((8)) with n = 3 and 



£2 As 



(1 + cos0)gii„ In 



47ro-eAB 



1 2 



(35) 



Here, 1 < ei < 3 and £2 ~ 80 are the dielectric constants of air (or oil) and water, respectively. As « 0.7 nm is the 
Bjerrum length in water at room temperature, k^j^ ~ 1 /^^i is the Debye screening length in water, e is the positive 
elementary charge, CTc is the charge density at the surfaces of the particles in contact with water, is the contact 
angle at the particle-interface contact line, and 0.1 ^ g\in{0,i^D) ^ 1 is a factor of geometrical origin. For strongly 
charged colloids one typically has CTc ^ 0.5 e/nm^, while water with a salt concentration above 10~^ M has a screening 
length kJ)^ below a few nanometers. Under these conditions, the factor [(1 +cos6').gii„ ln{4'Kac\B / ^kd)]^ is of the 
order of the unity, and in the following quantitative estimates we replace it by 1. Furthermore, we take £2/61 « 40 as 
appropriate for an oil-water interface. For these values of the parameters one has 



17.22 / R 



For this expression Fig. [8] shows Jeans' length 1/K (see Eq. 
values q of the mean interparticle separation in units of R. 



/im 



2/3 



(36) 



33)) ) as a function of the particle radius R for various 



3. Repulsion between induced electric dipoles 

In the experiment described in Rcf. [66j , the dipoles induced by the external electric field give rise to an interparticle 
repulsion described also by Eq. © with n = 3 and 



eoR^E^ 
37^ 



■(£2 + ei)ip, 



(37) 



where the factor ip depends on the dielectric constants and the height of the particle at the interface (see Fig. 9 in 
Ref. [l^l, where the factor -0 is called fd)- For the experiment in Ref. one has ei = 1, £2 ~ 2.87, and « —0.019 
so that 
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E 



10^5 V/r 



4/3 



R 
/im 



(38) 



at room temperature. For this expression. Fig. [HI shows Jeans' length 1/K (see Eq. ([55)) ) as a function of the particle 
radius R. 



(b) 

3 mm 




FIG. 8: (color online) Jeans' length for buoyancy induced monopoles (see Eq. (126^ *) and (a) for the ideal gas equation of state (see 
Eq. ((STJ) (also valid for a gas of hard disks, see the discussion after Eq. p4|l '). or (b) for the equation of state corresponding to 
dipolar repulsion between electrically charged particles estimated by taking charge renormalization into account (see Eqs. (|33p 
and (|36l) ). q is the mean interparticle separation in units of R (Eq. ((24])). The dashed, horizontal line indicates the value of 
3 mm corresponding to a typical value of the capillary length A. Capillary-induced clustering is possible below this line (see 
Fig. EJ. (The curve for g = 3 in the inset stops at _R « 4 fim due to the onset of freezing.) 
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FIG. 9: (color online) Jeans' length 1/K for buoyancy and electric field induced monopoles (see Eq. (|29[) ) and for the equation 
of state corresponding to electric field induced dipolar repulsion (see Eqs. (|33p and ([38])). In (a) Jeans' length is reported for 
various values E of the electric field in units of 10*^ V/m and for a fixed mean interparticle separation g = 30 (in units of R). 
The divergence of 1/K at i? sa 15.7 nm for E = 0.5 corresponds to a vanishing total capillary monopole. In (b) various values of 
q for a fixed electric field = 2 x 10'^ V/m are studied. The dashed, horizontal line indicates the value of 3 mm corresponding 
to a typical value of the capillary length A. Capillary-induced clustering is possible below this line (see Fig. [S)). (When a curve 
stops, it does so at the corresponding onset of freezing.) 



4- Repulsion between induced magnetic dipoles 



In the experimental setup with superparamagnetic particles of susceptibility Xnn an external magnetic field H 
induces a magnetic moment, which in a spherical particle of radius R is given by 

m = XmH. (39) 

In turn this leads to a dipolar repulsion described also by Eq. ([8|) with n = 3 and 

^ (40) 

In the experiments described in Refs. 0, [13 the susceptibility is Xm ~ 1-7 and the magnetic field H ranges typically 
between 10^ A/m and 10'^ A/m, which corresponds to 

at room temperature. (For comparison, Earth's magnetic field has a strength of about 40 A/m and sets a lower bound 
on the value of H achievable in the laboratory, unless the magnetic field is generated in a specific configuration so as 
to counterbalance Earth's field.) Figure [TU] shows Jeans' length 1/K (see Eq. ([55]) ) as a function of the particle radius 
R. (As remarked in Subsec flV A3[ the contribution to the capillary monopole due to the magnetic field is negligible.) 



V. DISCUSSION 



Our results show that there is a range of parameters for presently accessible experimental setups which allows one 
to observe the instability of colloidal monolayers at fluid interfaces driven by capillary attraction: particle size in the 
micrometer range, Jeans' length well below the capillary length, and Jeans' time ranging from minutes to weeks. This 
provides theoretical evidence for the possibility that capillary attraction can lead to relevant aggregation effects at 
submillimeter length scales in spite of its relative weakness at the mean distances between the particles. 

In view of the estimates derived in Sec. IIVI it is conceivable that many experiments carried out so far happen to 
operate in the range of parameters within which the homogeneous configuration is stable against capillary-induced 
clustering. For charge-stabilized colloids (see Fig. [8jb)) the particles employed experimentally are usually too small 
(not larger than a few micrometers). In the experiments carried out with superparamagnetic particles [l^ (see 
Fig. ITU)) , the particles arc too small (R w 2.5 /xm), the magnetic fields too large, and it is likely that the capillary 
monopole is also too small (because the particles are completely submerged in water, albeit close to the interface). In 
Ref . ^661] , which reports experiments in an electric field with spherical particles ranging in radius between 2 /xm and 
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FIG. 10: (color online) Jeans' length 1/K for buoyancy induced monopoles (see Eq. (|26p ) and for the equation of state 
corresponding to magnetic field induced dipolar repulsion (see Eqs. (|33|) and (|4ip ). In (a) various values of the reduced 
mag netic field H = H/ilQ^ K/m) are studied for a fixed mean interparticle separation g = 30 (in units of R). In (b) Jeans' 
length is shown for various values of q at the fixed field strength H = 40A/m corresponding to Earth's magnetic field. The 
dashed, horizontal line indicates the value of 3 mm corresponding to a typical value of the capillary length A. Capillary-induced 
clustering is possible below this line (see Fig. [31). (When a curve stops, it does so at the corresponding onset of freezing.) 

77 /im, there is a brief remark on the clustering of particles by capillary attraction in the absence of the electric field: 
In these experiments the particles were in close contact in the clustered states (with a mean interparticle separation, 
as inferred from the photographs, of g ^ 2 (in units of R)). corresponding to states in the solid phase. The lack of 
information on the pre-clustered mean particle density renders our results in the fluid phase of limited use for the 
comparison with the experimental observations. Nevertheless, assuming a hard-disk equation of state for the particles 
in the absence of the electric field, Fig.|8l[a) indicates the formation of a cluster for the range of particle sizes employed 
in this experiment, with collapse times spanning several orders of magnitude (see Fig. EUa)). The only quantitative 
experimental studies of the clustering instabilit y w e are aware of which do provide data amenable to comparison 
with our calculations can be found in Refs. [isl. |69||. According to the interpretation of those authors, temperature 
inhomogeneities at the interface set the particles in tomotion, the clustering of which is initially driven by capillary 
forces; further restructuring of the emerging formations inside the clusters involves short-ranged forces. This latter 
feature is beyond the scope of our model and we are not in a position to judge the influence of the temperature 
inhomogeneities. We simply note that, for buoyant particles of size R 35 ^m spread with an initial areal density 
corresponding to 9 « 10 [4^1, Fig. [8] predicts indeed capillary-induced clustering and Fig. |6l[a) yields Jeans' time of 
the order of a few minutes, in good agreement with the reported characteristic times (sec Fig. 3 of Ref. [i^). In 
summary, according to our estimates it seems to be possible to perform experiments within an appropriate range of 
controllable physical parameters which would promote the occurrence of the capillary-induced instability in a variety 
of conditions and which would allow the systematic study of its dynamical evolution. 

The tempting question arises whether the so far unexplained interparticle attraction and the ensuing micron-sized 
clusters we referred to in Sec. [J and which are reported by various groups can be understood within the physical 
picture we have presented: a cluster would be held together against repulsion and thermal diffusion by the collective 
capillary attraction and it would be the final, equilibrium state of a capillary-induced collapse. The answer is negative. 
First, there is a dynamical counterargument based on our theoretical finding that a capillary-induced clustering would 
require the almost simultaneous collapse of spontaneous density fluctuations of all sizes (see Eq. (|23p ): shortly after 
the formation of the actually observed micrometer-sized clusters, there should also arise many other clusters of larger 
sizes and all particles would eventually gather in a single, large cluster. Such a phenomenon has not been reported. 
Secondly, there is a static counterargument based on the qualitative reasoning put forward in Subsec. Ill Al which 
states that the equilibrium size of the cluster is of the order of Jeans' length associated with the cluster density: as 
follows from Fig. [51 there is no range of realistic values of the parameters for which Jeans' length is of the order of 
the observed cluster size, typically tens of micrometers at most, i.e., less than 10~^ mm. In other words, the observed 
clusters consist of too few particles in order to be able to build up a collective capillary attraction of relevant strength. 

One has also observed colloidal crystals spanning the whole system, which has a typical size L in the millimeter 
or centimeter range, i.e., comparable with the capillary length A or somewhat larger (see, e.g., Ref. [HI, [6^). One 
could try to explain this particle distribution as a very large cluster, self-confined by its own capillary attraction. 
Within our model, this would correspond to a clustered state characterized by a cluster Jeans' length Lj of the order 
of A (see Subsec. Ill Ap . Although our results pertain to the fiuid phase of the 2D colloid, one could nevertheless use 
them as rough estimates for the solid phase, given the relatively weak squared-root dependence of Jeans' length on 
the compressibility (sec Eq. P^ ). Thus, Figs. [51 and [51 provide evidence that the condition Lj ~ A could be easily 
satisfied. Actually, Ref. [6y| mentions briefly the interpretation of the occurrence of large clusters as being due to 
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capillary attraction. 

Recently, Pergamenshchik 4J] has correctly pointed out the enhancement of the pairwise capillary attraction due 
to a collective effect involving many particles. His work addresses only the equilibrium configuration of clusters of a 
size much larger than the capillary length and in the solid phase. His somewhat involved calculations can be put in 
our present context as follows: The solid phase is described effectively by the equation of state of a harmonic solid 
of soft spheres, which for the repulsive potential given by Eq. (jH]) corresponds to a pressure p{g) oc at high 

densities [g^- Integration of the equilibrium condition expressed by Eq. pil) yields g oc {U — [/q)^^": where Uq is an 
integration constant. Inserting this result into Eq. © one recovers Eq. (39) of Rcf. [i^l, which is his central result 
(notice that the exponent n in our notation (see Eq. ([8|)) corresponds to n — 1 in Pergamenshchik's notation (see 
Eq. (31) in Ref. These considerations show that the approach used in Ref. [i^l is contained in the theoretical 

framework presented here. 

However, in view of the discussion given above his claims that such a collective effect explains the clusters observed 
so far in experiments appear to be unsubstantiated. First, the analysis of an infinitely extended cluster misses the 
explicit size-dependence of the capillary energy (see Eq. (|4])), which makes his results unapplicable to small clusters. 
Secondly, his actual application of these results to experiments considers only either an ideal gas or a gas of hard 
disks. As we have shown in our analysis, these two models do not give rise to significant differences among them 
but both represent inadequate approximations for generic conditions in typical experimental systems. This difference 
can be understood in terms of the non-capillary contribution to the energy (see Eq. (O): an ideal gas or a gas of 
hard disks contributes kinetic energy only, which is of the order of the thermal energy T per particle (in units of the 
Boltzmann constant), whereas for the interparticle potential given by Eq. ([5]), the virial theorem yields an energy per 
particle eshort/T" = 1 + {2/n){p/{gT) — 1), which according to Fig. [TJa) can be at least « lOOT in the solid phase. 

The analogy of the model presented here with the 2D evolution of a self-gravitating fluid is not complete: the 
capillary attraction is, unlike gravity, screened beyond the capillary length, and the temporal evolution is ruled by an 
overdamped dynamics (see Eq. (|13p . amended in general with the effect of hydrodynamic interactions at sufficiently 
high densities) rather than by the inertial, Newtonian dynamics for gravitating particles. This poses the question as to 
which extent the gravitational phenomenology can be reproduced by colloids at a fluid interface. Our study provides a 
partial answer, in that it demonstrates the existence of a clustering instability which is analogous to Jeans' instability. 
But this is still far from being a complete and systematic comparison, which could provide the tempting picture of the 
feasibility to study "cosmology in a Petri dish". In this context, it would be interesting to investigate the equilibrium 
conflguration of clusters and their stability beyond the simple qualitative analysis we have presented in Subscc. Ill Al 
i.e., by solving Eqs. ^ and ([TT|) . Such a study would be complementary to the analysis of the dynamical aspects 
we have addressed here. The form of the capillary attraction also leads to a possible analogy with two-dimensional 
vortices, the similarity of which with a self-gravitating system is also well known. However, how deep and useful this 
latter analogy can be is still a matter of study (see, e.g., Ref. [t^ and references therein). 



VI. CONCLUSION 



We have presented a mean-fleld model for the evolution of the density of colloidal particles at a fluid interface 
driven by its own capillary attraction. In spite of the weakness of the capillary interaction at the mean distance 
between particles, its non-integrable character at submillimeter length scales enhances its effect on the evolution of 
collective modes. We have demonstrated that if the characteristic Jeans' length (see Eq. ([TO)) ) of a homogeneous 
distribution is sufficiently small (see Fig. [3]) the system can be unstable with respect to long-wavelength density 
perturbations under the action of capillary attraction (see Fig.|4]). Beyond this linear stability analysis, we have also 
solved the nonlinear evolution of radially symmetric density perturbations (see Fig. [TT]) in the so-called cold-collapse 
approximation (within which the dynamics is driven by capillary attraction only) (see Figs. [T^- I14p . which predicts 
a typical time of collapse of the order of Jeans' time (see Eqs. (|T^ and (|C13p ). By computing Jeans' length and time 
for presently accessible experimental setups we obtain clear predictions about the range of parameter values within 
which the instability could be observed. Jeans' time (see Fig. [5]) depends on the strength of the capillary monopole. 
To this end we have considered the monopole to be induced either by buoyancy or by an external electric field. Jeans' 
length depends additionally on the equation of state of the 2D gas via its compressibility. In this context we have 
studied an ideal gas (see Fig. [8ja)) and systems with a dipolar interparticle repulsion (see Eq. (|8]) for n = 3) induced 
either by electric charges on the particles (see Fig.[5|b)), an external electric field (see Fig.[5|), or an external magnetic 
field (sec Fig. llOp . In most experiments performed so far, the physical parameters lie in the region of stability, but 
they appear to be easily tunable into the instability regime. The relatively weak dependence of Jeans' length on the 
equation of state (via the square root of the compressibility, see Eq. ([TO|). renders the capfllary monopole to be a 
possibly more convenient parameter for tuning Jeans' length and time. Experiments with 2D colloids exposed to an 
external electric field seem particularly promising in this respect, because the field provides a simple way of controlling 
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the vertical pull on the particles generating in turn the mediating interfacial deformation. Similarly, a vertical force 
on superparamagnetic particles in an external magnetic field could be created and controlled via gradients of the 
magnetic field which pull on the induced magnetic moments. 
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Appendix A: Functional formulation 



The mathematical model defined by Eqs. ^ and can be reformulated in terms of the functional T[g, U] ~ 
J^cap + -^gas + -^intor, Consisting of three contributions. The first one is related to the capillary deformation in the 
small-deformation limit: 



cap •= ^7 / dA 



(Al) 



The second term is the free energy functional (within local approximations) of the 2D gas of particles, 

J-gas J dA fg3s(r, g), fg3s(T, g) = Tg[ln{A^g) - 1] + fe,(r, g), (A2) 

where the free energy density fgas{T, g) is the sum of the ideal gas contribution (A is de Broglie's thermal length) and 
the excess free energy fex due to the repulsive short-ranged forces. Finally, the third term represents the interaction 
between the particles and the interfacial deformation: 



/ / dAgU. 



'Winter • 

The mean-field equation ^ for the interfacial deformation follows from the extremal condition 

5U{v) ' 

while the diffusion equation (|14p can be expressed in terms of a relaxation-type dynamics. 



dg 



gV 



_5T_ 



upon using the thermodynamical identity 



pie) 



2 d 
dg 



fgasC^. g) 



g 



at constant T. 



(A3) 



(A4) 



(A5) 



(A6) 



In principle, F can be viewed as an effective functional reformulation of the problem, although it can be associated 
with an actual free-energy functional for the physical system. In this case, the mean-field approximation enters 
via the simplified form of Eq. (|A3p : a complete description of the particle-interface interaction should take into 
account the finite size of the particles (rather than a point capillary monopole), their shape, the corresponding surface 
energies (see, e.g., the free-energy functional introduced in Ref. [38i|). and a thermal noise contribution to Eq. (jASp . A 
procedure for incorporating the hydrodynamic interactions into the functional formulation has been proposed recently 
in Ref. 
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Appendix B: Equation of state from numerical simulations 



We investigate a two-dimensional fluid governed by the pair potential given by the power law in Eq. ([8]) with n = 3. 
The equation of state of such a fluid depends only on the dimcnsionlcss parameter C,"^ p [s^]- Thus it is advantageous 
to introduce the new length scale a = '^j ^fp such that 

^-(cvf/^(5)^ (Bi) 

In the simulation one fixes b as the unit of length and varies the prefactor [C^ pf"^"^ . The pressure can be determined 
from the virial expression 

1 



i<j 



P^pT+—0 . f (r« - r,) • (r. - r,)) (B2) 



where A is the area of the simulation box and f = — Vi'rop is the intcrparticle force. In our case 



f(ri - Tj) • (rj -Vj) = 3 V'wi.opCi'j - r^) , (B3) 



/ ^ --JJ V-4 •■31 ~ " / ^ 

i<j i<j 

SO that the pressure is related to the excess internal energy If^^ of the fluid according to 

O JTCX 

The simulations have been carried out with N = 400 particles in a square simulation box with reduced side length 
L ~ 20 (in units of a). Periodic boundary conditions were applied. The internal energy of a given configuration of 
the N particles was obtained by summing over {21 + 1) x {21 + 1) boxes (i.e., the central simulation box plus image 
boxes around it) because the pair potential decays slowly (we chose 1 = 3). U^^ was determined as an average over 
10^ sweeps (one sweep corresponds to attempted moves for every particle in the box). The instantaneous internal 
energy was determined for each sweep once. 



Appendix C: Radially symmetric cold collapse 




FIG. 11: Evolution of a radially symmetric enhanced density. An infinitesimally thin ring of particles of initial radius Rq 
collapses to a ring of radius R{t) < i?o at a time t > 0. By construction the ring follows the motion of its particles and it can 
be identified by its Lagrangian radial coordinate Rq (and R{t) is the Eulerian radial coordinate of the ring). Therefore, the 
number of particles inside the disk encircled by the ring is constant if different rings do not cross (one cannot exclude a priori 
that a ring is overtaken by another one because different rings may collapse with different speeds). 



1. General dynamics of radially symmetric cold collapse 



A homogeneous configuration is characterized by a density g/^ and a mean-field interfacial deformation = 
jy? QhH given by Eq. We study the evolution of a radially symmetric density perturbation so that 

Q{v,t) = Qn + ^e{r,t). (CI) 
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The perturbation SU of the interfacial deformation is given by the solution of Eq. (|T6)) . (Note that this latter equation 
is not a linearized approximation requiring Sg and 6U to be small, because Eq. is linear; the nonlinearity of the 
dynamics enters via Eq. ([M]) .) We consider a sufficiently localized density perturbation, i.e., Sg vanishes sufficiently 
fast at infinity, so that the solution in the limit of large capillary length A follows immediately from the gravitational 
analogy: the capillary force on a particle is (compare Eq. ([3])) 



fV(dU) = e,., 

2717?' 



(A oo) 



(C2) 



where 



dN{r,t) := 27r / dss5g{s,t) 



(C3) 



is the excess number of particles in a disk of radius r concentric with the perturbation. The density field follows 
from solving Eqs. (|12p -(|14 p . which in the cold collapse approximation, i.e., after neglecting any force other than the 
capillary attraction, reduce to (see Eq. (|C2[) ) 



d{Sg) 
dt 



-V • (vg), 



(C4) 



v = r/v(5c/). 



(C5) 



An exact solution is possible by resorting to a description of the flow in terms of Lagrangian coordinates. Consider 
the thin ring of radius R(t; Rq) formed by the particles which were initially {t = 0) at a distance Rq from the center, 
i.e., the ring is defined through its Lagrangian radial coordinate Ro (see Fig. [TT|) . Equation ()C4p implies that the 
number of particles encircled by the ring is conserved^: 



rR i-Ro 

N{R,t)-~ 271 dss g{s,t) ^ No{Ro) ^2tt ds s go{s) 
Jo Jo 



(C6) 



where 7Vo(i?o) is the initial number of particles encircled by the ring of radius i?o- (Since, at fixed t, N{R,t) is a 
monotonous function of R, Eq. (jC6p provides a relation Ro = Roit; R) which can be inverted, R = R{t; Ro), with the 
geometrical meaning "radius of the ring at time t which had a radius Rq at time t = 0".) This expression allows an 
explicit computation of the density field. 



N I dN 
'^''^'^=2^ OR 



1 dNo dRo 



2T:RdRo dR 

The evolution of the ring radius i?(t; Rq) is determined by the flow velocity, 



1^^°^^°) -dR 



dR 
dt 



e,. V, 



(C7) 



(C8) 



Ro 



with V given by Eq. (jC5|) . By inserting the solution (jC2|) and making use of particle conservation, 6N{R,t) = 
N{R, t) — TiR^gh = No{Ro) — T^R^gh (see Eq. (jC6P '). one finally arrives at an ordinary, first-order differential equation 
for R{t]Ro): 



dR 
-dt 



Rq 



2'K-iR 



[No{Ro)-^R^gh] 



(C9) 



In terms of the characteristic time T — l/i^Pgh) of the homogeneous background configuration (sec Eq. (fTn|) ) one 
can write 



dt \Ro 



R 
Ro 



Qh' 



(CIO) 



In this appendix, a subindcx denotes evaluation at the initial time t = 0. 
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where 

MR.):^^ (Cii) 

is the average initial density of the disk of radius Rq. The solution of this equation with the boundary condition 
R{t = 0; i?o) = Ro reads 



i?(t;i?o) = i?oi/-+ (l--)e*/r. (C12) 

^ eh V ShJ 

Figure [1^ plots the evolution R{t) for the two possible, qualitatively different cases: If Qh > Qq the ring expands. 
(When R reaches A the model ceases to be valid because then one must take into account the screening of the capillary 
interaction neglected in Eq. (jC2|) .) In the other case, Qh < go, the ring contracts and i? = is reached at a time 

Tcoii = -Tin (l - ■ (C13) 

Obviously, true collapse is prevented by a finite compressibility which has been neglected in the above calculations (see 
Eq. (jCSp ) and a spatially extended equilibrium configuration will emerge. However, the unambiguos identification of 
a ring by its Lagrangian radius Rq would no longer hold if two different rings would coincide before the collapse at the 
center (a phenomenon customarily called shell crossing in 3D in the cosmological literature). If this would happen, a 
density singularity would appear, maybe regularized by the repulsive forces present. It can be shown that a sufficient 
and necessary condition to avoid ring crossing is that the time of collapse 7^oii(-Ro) of the ring characterized by its 
initial radius Rq grows with this radius. In view of Eq. (jClSp . this is equivalent to the condition 

i.e., the initial perturbation inside a disk looks more and more rarefied as the disk radius is taken larger and larger. 




t/T 

FIG. 12: (color online) Evolution in time of the radius R{t;Ro) of a circular ring during cold collapse (Eq. (|C12[l '). Two 
qualitatively different cases (expansion and collapse) can be distinguished depending on the initial average density in the 
area encircled by the ring compared to the background density Qh- The curves correspond to Qo/ Qh = 0.5 for expansion and 
Qo/Oh = 1.5 for collapse. 



2. Example of cold collapse for a steplike overdensity 

As a simple application of these results we consider a top-hat enhanced density, described by the initial profile 

Qoir) = + 5go e(Lo - = ( ^'^ t \ (Sqq > 0), (C15) 
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where 8(x) is Heavyside's step function. We assume that the initial radius fulfills K ^ Lo ^ A, so that the 
simplifying assumptions of the previous calculations hold. For this profile, one has (Eq. (|Clip ) 



Qh + Sgo, 



QoiRo) 



Qh 



Ro 



(C16) 



^90, Lq < Rq, 



SO that Qh < go and, according to Fig. [T^l there is a collapse for all Lq > 0. The time of collapse (sec Eq. (|C13P ') 

2 



rcoii(i?o) = rhJ 1 + 1^ 



(C17) 



is plotted in Fig. 1131 the whole interior of the overdensity (Rq < Lq) collapses simultaneously, i.e., T^oii is independent 
of Rq- The evolution of the density profile g{r,t) follows from Eqs. (|C7[) and (jC15[) : by inserting Eq. (|C16[) into 
Eq. (jC12[) one obtains upon differentiation 



R OR 



_ j (i?/i?o)^, ^0 < Lq 
1, Lq < Rq 



Rq dRo 

so that the density profile preserves the top-hat shape during the evolution (see Fig. [T4|): 



(C18) 



g{r,t) = 



Lo 



[gh + Sgo), r < L{t) 
L{t) < r 



(C19) 



gh 



where 



Lit) = Loi 



gh + 5go 
gh 



[1 



.(t-TooiO/rl 



(C20) 



is the time-dependent radius of the region with enhanced density, i.e., Eq. (|C12p evaluated for the particular case 
Rq = Lq by using Eqs. (jClSp and (|C16[) . Consistently with this geometrical meaning, L{t) vanishes at the time of 
collapse t = 7^ou(-Ro = ^o) of this region, as given by Eq. (|C17|) . (Upon deriving Eq. (|C19|) one notices that the 
region with enhanced density can be characterized by the inequality Rq < Lq or equivalently by r < L(t), since 
r = R{t; Rq).) The sharp jump in the density profile is of course an idealization in the limit of infinite compressibility. 
The jump would be actually smoothed out by the pressure forces, but on a length scale much smaller than the spatial 
extension of the density enhancement. 

As the perturbation is compressed, in the 2D fluid the pressure increases until it is able to counterbalance the total 
capillary force and stop the collapse. One can derive a simple estimate for the smallest size Lmin of the compressed 
cluster as follows: at the maximum density £>max ^ gh-, corresponding to a cluster radius imin <^ Lq <^ X, a total 
pressure force^ of the order of Fp ~ 27rLminP(gmax) opposes the total compressing capillary force, which is of the order 
of Fcap ^ (/7V)^/(27r7Lniin)- (The cluster consists oi ^ N particles and the capillary force on each of them is given 



by Eq. ^ with 6N 
relationship 



Qh) ~ N due to Pmax ^ gh-) Therefore, in equilibrium Fp ~ Fcap, giving the 



Pigmux) 



mm t/maxy 
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(C21) 



Upon introducing Jeans' length i^^iustcr °f cluster, i.e., Eq. evaluated at the cluster average density g),, 
this can be rewritten as 



luster-' 



0^ 4p(glmax)K(£'max)- 



(C22) 



^ Note that in 2D pressure is a force per unit length. 
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FIG. 13: (color online) Dependence of the time of collapse on the ring radius for a top-hat enhanced density distribution. The 
plot corresponds to the case 5 go — O.bQh (see Eq. (|C17|l 'l. The time scale is set by Jeans' time T (Eq. (US}). 
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FIG. 14: (color online) Time evolution of the density profile of an initial top-hat enhanced density distribution (Eq. (|C19|) V 
The plot corresponds to Sqo — 0.5gh- The time scale is set by Jeans' time T (Eq. (|19p ). 



For polytropic equations of state, p{g) oc g^, or more generally for a simple fluid far from any phase transition, the 
right hand side of this expression is typically of the order of unity, so that the final equilibrium size of the cluster 
will be comparable to its Jeans' length. This rough estimate agrees with the results obtained from the qualitative 
reasoning given in Subsec. Ill Al 
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